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1 Introduction 



From time to time the detailed study of the results of the analytical multiloop calculations 
allow one to reveal the existence of the internal symmetries and of the definite properties of 
the gauge theories under investigation. In fact these typical features of the concrete models 
can be hidden in the explicit expressions of the coefficients of the perturbative series for the 
renormalization group quantities. The classical example of the immediate influence of the 
outcomes of the analytical calculations to the further development of the understanding of the 
structure of the perturbative series as the whole is provided by the evaluation of the 3-loop 
contribution to the /^-function of the N = 4 supersymmetric Yang-Mills theory [1]]. Indeed, 
in Ref. [1] it was found by the direct diagram-by-diagram calculation that the corresponding 
/3-function is zero at the 3-loop order. In its turn, this foundation pushed ahead theoretical 
works, which resulted in the formulation of the proof of the validity of this interesting property 
of the N = 4 supersymmetric gauge model in all orders of the perturbation theory (see e.g. 
Ref. [2]). Therefore, in spite of the fact that the calculation of the related /3-function coefficients 
presumes the introduction of the regularization (definitely speaking, supersymmetric one) and 
thus the renormalization procedure, the property of the conformal symmetry is preserved in 
this model. 

However, it is known, that in such realistic gauge theories, as massless QCD or QED, the 
procedure of the renormalization is disrupting the initial conformal invariance and gives rise 
to an anomaly in the trace of the energy- momentum tensor || . Its explicit expression M 
shows that the factor j3{a)/a is the measure of the breaking of the conformal invariance within 
the framework of the perturbation theory expansion in the coupling constant a, where in QCD 
and QED we will normalize a as a s = a s /n and a = a/n respectively. This means, that in 
QCD and QED the conformal symmetry can be effectively restored only in the vicinity of the 
hypothetical perturbative fixed points, which satisfy the condition (3{a*) = 0. In view of this 
it is rather interesting to understand whether there are any manifestations of the properties 
of the initial conformal symmetry and its violation by the procedure of the renormalization in 
the structures of any perturbative series. 

By what discussed beyond we will try to convince the readers that the answers to these 
questions are positive. We will show that the made in Ref. [5] careful analysis of the analytical 
structure of the perturbative QCD predictions to the certain characteristics of the e + e~ — > 
hadrons and deep-inelastic scattering processes allowed to reveal the existence of the definite 
already proved |J and still non-proved relations between the coefficients of the perturbative 
series for the Adler D-function of the non-singlet axial currents (or vector currents) and 
the polarized Bjorken sum rule (or the Gross-Llewellyn Smith sum rule). The first quantity 
was evaluated analytically at the next-to-leading order (NLO) of perturbative QCD in Ref. [7] 
(for the identical result of the independent numerical and analytical calculations see Ref. [8] 
and Ref. [9] respectively) and at the next-to-next-to-leading order (NNLO) in Ref. [10] (for 
the identical result of the semi-independent calculations see Ref. [11]). The NLO coefficient 
of the perturbative series for the deep-inelastic scattering sum rules we will be interested in 
is known from the results of the analytical calculation of Ref. [12] (later on independently 
confirmed in Ref. [13]), while the NNLO corrections to the polarized Bjorken sum rule and to 
the Gross-Llewellyn Smith sum rule were obtained in Ref. [14] . 

Starting from the beginning of the INR multiloop analytical calculating project, launched 
in 1978-1979 0], nobody was expecting that there are any relations between the character- 



istics of the annihilation and deep-inelastic processes. However, in 1990, when it was already 
understood by the second and the third authors of Ref. [10], that the published results of the an- 



alytical calculations of the 4-loop contributions to the QED /5-function in the MS'-scheme [16 



and to the Adler .D-function in QCD |l7j are wrong (for the discussion see e.g. Ref.[18]), but 
the 4-loop QED results of Ref.[19] and the NNLO QCD results of Refs. [10,11,14] were not yet 
obtained, we were informed in the quite non-formal form p0| about the existence of the funda- 
mental Crewther relation |3|]. This relation is connecting in the conformal- invariant limit the 
anomalous 3-point function of the axial-vector-vector non-singlet currents with the product 
of the quark-parton expressions for the polarized Bjorken sum rule and the e + e~-annihilation 
Adler .D-function. However, in order to answer to the constructive question: "What is the sta- 
tus of the Crewther relation in QCD?" [^] it is necessary to go beyond the framework of the 
conformal-invariant limit and to analyze the structure of the perturbative QCD corrections to 
both sides of the Crewther relation. 

Here we will summarize our present understanding of the current status of the answers to 
this question using the results of the works of Refs. [5,6,21] and of the old work of Ref. [22], 
which became known to us only recently. 



2 What is the Crewther relation? 

Before starting the presentation of the QCD foundations of Ref. [5] let us briefly discuss the 
definite steps of the derivation of the Crewther relation ||, repeated also in Ref. [22]. In 
order to complete the analysis of Refs. [3,22], performed in the x-space, we will follow the 
studies of Ref. [6] and use the language of the momentum space. To our mind, this will 
allow to demonstrate more obviously some basic points, which were not previously clarified in 
Refs. [3,22]. 

Consider the 3-point function 

T^ip, q) = ij < 0\TA«(y)V^x)Vp C (0)\0 > e^ + ^dxdy = d abc T^(p, q) (1) 

where A a ^(x) = ?/>7 At 7 5 (A a /2)?/>, V®(x) = ipj^X 0, /2)ip are the axial and vector non-singlet quark 
currents. The r.h.s. of Eq.(l) can be expanded in a basis of 3 independent tensor structures 
under the condition (pq) = as 

+ i2(p 2 ,q 2 )(q a e^p T p p q T - q^ mpT p p q T ) (2) 

+ iz(p 2 ,q 2 )(p a t^PTP p q T + PP e mprP p q T ) ■ 

Taking now the divergency of axial current one can get the following relation for the invariant 
amplitude £,i(q 2 ,p 2 )'- 

qpT^ aP (p, q) = e liapr q p p T ^ 1 (q 2 J p 2 ) (3) 
while the property of the conservation of the vector currents implies that 

/iim|p2_ 00 p 2 ^3(g 2 ,p 2 ) = -ii(q 2 ,p 2 ) (4) 

(see Ref.[23] for the discussions of the details of the derivation of Eqs.(2)-(4)). 

In order to clarify the meaning of the second invariant amplitude, namely £,2{q 2 ,P 2 ), let 
us first define the characteristics of the deep-inelastic processes, namely the polarized Bjorken 
sum rule 

Bjp(Q 2 ) = l\g7(x,Q 2 )-gr(x,Q 2 )]dx = h^\C Bjp (a s ) (5) 
Jo 6 gv 

and the Gross-Llewellyn Smith sum rule 

GLS(Q 2 ) = \ I' F» p+T7p (x, Q 2 )dx = 3C GLS (a s ) . (6) 
I Jo 



The coefficient function Csj P (a s ) can be found from the operator-product expansion of two 
non-singlet vector currents JTT 



where 



i J TV^VjftSje^dx^^ « C^ a p C A c p (0) + other structures (7) 

C a0p C ~ id abc e a pp a j^C B jp{a s ) . (8) 

and P 2 = — p 2 . In the case of the definition of the coefficient function of the Gross-Llewellyn 
Smith sum rule one should consider the operator-product expansion of the axial and vector 
non-singlet currents 

i J TA^V^Oy^dxl^^ « C%£v a (Q) + other structures (9) 

where 

CJS ~ ti ab ^^C G Ls{a s ) (10) 

and Q 2 = — g 2 . The third important quantity, which will enter into our analysis, is the QCD 
coefficient function C^ s (a s ) of the Adler P-function of the non-singlet axial currents 

D NS (a s ) = -12vrV ^Nstf) ~ C» a {a.) (11) 

where Uns (q 2 ) is defined as 

i J < 0\TAl(x)Al{0)\0 > e^ x dx = 5 a \g, u q 2 - q^)H N s{q 2 ) • (12) 

At this point we will stop with definitions of the basic quantities and return to the consid- 
eration of the 3-point function of Eq.(l). Following Ref. [3] one can apply to this correlation 
function an operator-product expansion in the limit \p 2 \ » |g 2 |, p 2 — > oo, namely expand 
first the T-product of two non-singlet vector currents via Eq.(7) and then take the vacuum 
expectation value of the T-product of two remaining non-singlet axial currents defined through 
Eq.(12). It was shown in Ref. [6] that these studies imply that 

6(<? 2 ,P 2 )l| P 2 Hoo -> ^C B jp(a s )Yl NS (a s ) (13) 

and thus 

q 2 ^M^^%^oo^j 2 C BjP ^s)C^ s {a s ) . (14) 

Equations (13), (14) reflect the physical meaning of the invariant amplitude ^2(? 2 ,p 2 ) and 
should be considered together with the relations for the invariant amplitudes £,i{q 2 ,p 2 ) (see 
Eq.(3)) and £ 3 (? 2 ,P 2 ) (see Eq.(4)). 

On the other hand, it was shown in Ref. [24] that in a conformal invariant (c-i) limit the 
three-index tensor of Eq.(l) is proportional to the fermion triangle one-loop graph, constructed 
from the massless fermions, namely that 

T%>(P, <l)\-i = d abc K(a s )Al- l °° p (p, q) . (15) 



A 



In other words, in a conformal invariant limit one has 

(T{M = K(a s )er l °° P (q 2 ,P 2 ) , 

£r(<rV) = K{a s )e 2 - l °° P {q 2 y) , (16) 

er(? 2 ,P 2 ) = K(a s )e 3 - loop (q 2 , P 2 ) • 

Moreover, in view of the Adler-Bardeen theorem p5f , which is insuring the invariant amplitude 
£i(q 2 ', p 2 ) , related to the divergency of axial current (see Eq.(3)), from the renormalizability, 
one has K(a s ) = 1. The 3-loop light-by-light-type scattering graphs, which were calculated 
in Ref.[26] and analyzed in Ref.[27], do not affect this conclusion. Indeed, in the case of the 
3-point function of the non-singlet axial-vector-vector currents they are contributing to the 
higher order QED corrections, while the QCD corrections of the similar origin are appearing 
only in the 3-point function with the singlet axial current in one of the vertexes, which will 
be not discussed here. 

In the case of the consideration of the 3-point function of the non-singlet axial- vector- vector 
currents the property K(a s ) = 1 is allowing to derive the fundamental Crewther relation 

C Bjp (a s (Q 2 ))CS s MQ 2 ))\c-i = l , (17) 

which should be valid in the conformal invariant limit in all orders of perturbation theory. The 
similar relation is also true for the coefficient function Cgls^s), defined by Eqs.(6),(9),(10) [p2[ . 
Indeed, considering first the operator-product expansion of the axial and vector non-singlet 
currents (see Eq.(9),(10)) in the 3-point function of Eq.(l), taking the T-product of the re- 
maining vector currents and repeating the discussed above analysis, one can find that in the 
conformal invariant limit the following identity takes place: 

C GL s(a s (Q 2 ))CUa s (Q 2 ))\c-i = 1 (18) 
where C^(a s ) is the coefficient function of the Adler Z)-function of two vector currents. 



3 The QCD generalization of the Crewther relation 

It is well known, that the calculations of the perturbative theory corrections of the Green 
functions in the renormalizable quantum field models face the necessity of the introduction of 
the concrete regularization of the ultraviolet divergencies and of the subsequent application 
of the renormalization procedures, which as already mentioned, are breaking the conformal 
symmetry of the massless free theories. The immediate consequence of the application of this 
modern perturbative machinery is the appearance of the renormalization group /3-functions, 
which are governing the energy behavior of the running coupling constants and are responsible 
for getting out from the scale-invariant limit. In QCD the /^-function 

Q 2 -^ = /3(a,) = -£/K +2 (19) 



was analytically calculated at the 3-loop order in Ref.[27] and confirmed in Ref.[28] in the 



framework of the dimensional regularization [30] and the class of minimal subtractions schemes [31 



In our further discussions we will be interested in the expressions of the first two renormalization- 
scheme invariant coefficients /3q and /3i, which are expressed through the Casimir operators 
Ca, TfNf and Cp as 



A) = (f<?A-iW)± 



1 ' x (20) 

A = (%a%-fC A T f N f -4C F T f N f )± . 



E 



In general, it is possible to understand that 



• the perturbative expression for the QCD /3-function does not contain the terms propor- 
tional to Cp af (K > 1) and 

• the deviation from the conformal invariant limit (3{a s ) = is related to the fact that 
T f N f ^ and C A ^ 0. 

Moreover, in perturbative series for physical quantities the latter Casimir operators can appear 
only starting from the NLO. Therefore, to study the theoretical consequences of the property 
of the conformal symmetry breaking in the massless gauge models it is necessary to consider 
the higher-order perturbative theory approximations for the physical quantities. 

In order to get the non-conformal variants of the Crewther relations of Eqs.(16),(17) it 
is necessary to consider the NNLO approximations of the corresponding basic quantities. 
Due to the existence of the dimensional regularization []3C[1 , class of the minimal subtractions 
schemes [BTj and the developments in the field of the creation of the multiloop calculat- 



ing methods [15], j32|, it became possible to obtain the concrete NNLO results for the Adler 



function of the electromagnetic quark currents [ID, [TT| by means of the classical symbolic 



manipulations program SCHOONSCHIP p3[ and for the QCD coefficient functions of the po- 



larized Bjorken sum rule and the Gross-Llewellyn Smith sum rule |14[ with the help of its more 
computer-educated younger follower FORM ||34|| . Another problem, solved in the process of 
the NLO calculations of CBj P (a s ) and CGLs( a s) of Ref.[12] and of the NNLO ones of Ref.[14] 
is the proper definition of the axial non-singlet current. 

Indeed, it is known, that within dimensional regularization the analog of the 75-matrix and 
therefore the axial Ward identities are not well defined. The most straightforward way, which 
is allowing to restore the axial Ward identities presumes the application of the additional finite 
renormalization of the axial currents P5| , which in the non-singlet case has the following 
form 

Al^>Z» 3 Al = Al . (21) 

In the MS-scheme the expression for Z§ s was calculated at the NLO level in Ref. [12] and 
NNLO level in Ref. [14]. It is possible to show, that this additional finite renormalization 
allows one to restore the validity of the Adler-Bardeen theorem in QCD for the anomalous 
triangle diagram of the axail-vector-vector non-singlet currents (for the related discussions 
see Refs. [36-38]). Moreover, using Eq.(21) one can find, that the NNLO approximation of 
the QCD coefficient function of the Adler D-function of the electromagnetic quark currents 
C^) M {a s ) differs from the NNLO approximation of C^ s (a s ) only by the 4-loop light-by- light 
type diagrams, which have the singlet nature: 

CE M "(a.) = C% s '(a.) + Cf (a s ) . (22) 



The last scheme-independent contribution has the following form WQ, [TT 



where n c = 3 is the number of colours. In the MS'-scheme the NNLO analytical expression of 
Cp S (a s ) is known from the results of calculations of Refs. [10,11]: 



C% s (a s ) = 1 + hj F a s 



,123 11 ^ \ /ll 



S2° F + {^2~ T C(3) ) CfCa+ \-- + ^(3) )C F T f N, 



a 2 . 



+ 



+ 



^90445 



3456 



2737 
144 



69 
128 

C(3) 



C F +[ - 



55 
24 ( 



29 19 



(24) 



C(5) )C F C A + - - + -C(3) - 5C(5) )C F T f N 



485 + ^C(3) + §C(5) )C F C A T f N f + 



27 



9 



64 4 
/151 19 



C(3) )C F TjN 



The similar MS'-scheme result for the coefficient function CBj P (a s ) was obtained in Ref.[14] 
and reads 



+ 



CBjp{.CL 

3 



1 - -C F a s + 



21 23 1 

—Cl C F C A + -C F T f N f 

32 F 16 2 ' ' 



128 



C 3 P 



/1241 11 



V 576 



^((3))C 2 F C A + 



Z3535 3 Wo . 



6 



5437 + 55 ^V r 2 
W + 24 C(5) J CfCa 

C{$))c F C A T f N f -^C F T}N} 



(25) 



while the NNLO coefficient function CcLs{ a s) recieves the additional NNLO singlet-type 
contribution [14| 

C GLS (a s ) = C Bjp (a s ) + C^ LS (a s ) (26) 



where 



n 



1 192 8 SV ; ' 



(27) 



From the first superface glance to the pairs of the results of Eqs.(24),(25) (or to the related 
ones of Eqs.(22),(26)) one can think that there is nothing in common between the beautiful, 
but rather complicated expressions for the NNLO approximations of the characteristics of 
the annihilation processes (see Eqs.(22)-(24)) and the ones of the deep- inelastic processes 
(see Eqs.(25)-(27)). However, after constructing the NNLO QCD variants of the Crewther 
formulae of Eqs.(17),(18) and looking to these products more carefully the following relation 
was found |J: 



1 + 



C B]P {a s {Q 2 ))C* s {a s {Q 2 )) 



a 



S 1 C F a s + ( S 2 T f N f + S 3 C A + S A C F )C F a 



.2 



+ 0(a 4 s ) (28) 



where f3^ 2 \a s ) = —/3oa 2 — (3\a? s is the 2-loop expression for the QCD /3-function with j3 and 
Pi defined in Eq.(20) and Si, S2, 5*3 and 5*4 are the analytical numbers, which contain the 
transcendental Riemann (^-functions: 

Si = -f + 3C(3) , 

S 3 = -Jf + ^C(3) , (29) 
S * = W+ T C(3) - 15C(5) • 

The transformation to the non-abelian case of QED can be made by taking C F = 1, = 
and TjiVj = iV where A" is the number of massless leptons with identical charges. 

The authors of Ref. [5] noticed in Eqs.(28),(29) the following "seven wonders" of the gen- 
eralized Crewther relation, or to be more precise, of the NNLO non-conformal descrepancy 
A„_ c = (Cp S (a s (Q 2 ))CBj P (a s (Q 2 )) — 1) from the conformal-invariant relations of Eqs. (17), (18): 

1. The leading order terms cancels in A n _ c , which thus do not contain any CVa s -corrections. 

2. The NLO corrections give no Cpa 2 terms in A n _ c . 

3. The NNLO corrections give no C F a 3 s terms in A n _ c . These three foundations have led to 
the observation that in the so-called quenched limit the zero-fermion-loop abelian terms in the 
NLO and NNLO approximations of the polarized Bjorken sum rule and the Gross-Llewellyn 
Smith sum rule ( explicitely calculated in Ref. [12] and Ref. [14] respectivley) can be obtained 
by inversing the expression of the one-fermion-loop QED contribution to (a), which is 
related to the scheme-independent Baker- Johnson QED Fi-function, calculated at the 3-loop 
level in Ref. [39] (this result was confirmed in Ref. [40]) and at the 4-loop level in Ref. [19]: 

t~i / \ i i 3 21 2 31 g 

G Bjp{a>)\ quenched QED = 1 _ ~ a + — a — 77^77° 

(30) 



3 a + 4 fl2 32 a3 128 a4 



where the factor in the numerator is defined by the 1-loop contribution to the Fi-function. 

4. The NNLO light-by-light terms of Eqs. (23), (27) are cancelling in A n _ c (taking equal 
quark charges in Eq.(23) to obtain the contribution to C^(a s )). 

5. The NLO corrections give C f Cacl 2 and C F TfNf terms in A n _ c that are the same ratio 
as the Ca and TfNf terms in (3q. This reveals the exact factorization of the /^oa^-contribution 
in the NLO correction to A n _ c . 

6. The NNLO corrections to A n _ c can be expressed as the sum of the terms, proportional 
to (3qO? s and (3ia? s without introduction of the /^ag-terms, which can be expected from the 
simple power-counting arguments. 

7. At the NNLO level the /3ictf contribution to A n _ c occurs with the same coefficient that 
multiplies (3 a 2 term at the NLO. 

Before proceeding our discussions let us make several comments, related to the first three 
"wonders" mentioned above. It can be shown that these three properties are the conse- 
quence of the Adler-Bardeen theorem and of the discussed in Sec. 2 initial Crewther relation [0, 
which is valid in the conformal invariant limit TfNf = 0, Ca = 0. In fact they were already 



effectively discovered in Ref. [22], where the validity of the 3-loop analog of Eq.(30) was con- 
jectured. The result of the NLO calculations of CBj P {a s ) (see Ref. [12]), taken in the limit 
Cp = 1, TfNf = 0, Ca = 0, is confirming the prediction of Ref. [22] for the abelian zero- 
fermion-loop NLO contribution to this sum rule, while the results of the NNLO calculations 
of the F-function [19] and CBj P {a s ) 01 are extending the obtained in Ref. [22] 3-loop relation 
to the 4-loop level. Moreover, it is also possible to understand |J, that the properties of the 
cancellation of the Cj^af -terms (1 < K < 3) in the discovered in Ref. [5] NNLO analog of 
the Crewther relation (see Eqs.(28),(29)) can be generalized to the arbitrary order K of the 
perturbation theory. 

It is worth reminding here the physical meaning of the Baker- Johnson QED i<\-function. 
In the process of the study of the finite QED program it was proved that if the condition 
Fi(a*) = takes place, than the property of the finiteness of QED can be realized in the vicinity 
of the point a* |4"2]] . Thus, the explicite calculations of the coefficients of the perturbative 
series for the -F\-function can be considered as the important "experimental" ingredients of 
the analysis of the question of the existence (or, to be more precise, non-existence) of the 
perturbative ultaviolet fixed point of the massless QED. It is interesting to mention, that the 
3-loop analytical calculations of Ref. [39] revealed the cancellation of the £(3)-functions, which 
do appear at the intermediate stages of these calculations. As was mentioned later on in 
Ref. [40] without any reference, the fact that the 3-loop coefficient of the .Fl-function turned 
out to be rational can be related to the property of the conformal symmetry of the massless 
QED in the quenched approximation. 

It was also pointed out in the process of the discusions of the 4-loop results of the 1987-year 
calculation of the Fi-function fl6| , that it is rather doubtful, that the result of Ref. [16] (which 
was found later on to be in error) contained the £(5)-term El. Now we know, that in spite 



of the appearence of the C(3), C(4) and £(5)-terms at the intermediate stages of calculations 
of the 4-loop corrections to the -F\-function, they are indeed cancelling out in the ultimate 
correct result of Ref. [19]. Moreover, due to the appearence of the personal understanding of 
the place of the conformal symmetry in the derivation of the Crewther relation || 
we can agree with the comment of Ref. [40], that the perturbative theory expression for the 
-F\-function is constrained by the conformal invariant limit of the massless QED. We can also 
generalize the proposed in Ref. [40] hypothesis to arbitrary orders of perturbation theory and 
make the conjecture, that it is rather natural to expect that all coefficients of the Fi-function 
can be rational numbers and will not contain Riemann (^-functions. The arguments in favour 
of the yet non-proved in detail relations between the appearence of the Riemann ^-functions 
in other results of multiloop QED calculations, say the 2-loop and 3-loop analytical results 
for the anomalous magnetic moment of electron a e (see Ref. [44] and Ref. [45] correspondingly) 
and the theory of knots are given in the works of Ref. [46] . 

Let us now return to the discussions of the other four "wonders" of Eq.(28), which were 
discovered in Ref. [5]. In fact, since it is possible to understand (see Ref. [22] and the discussions 
in Sec. 2) that 

C B]P {a s {Q 2 ))C^ s {a s {Q 2 )) = C GLS (a s (Q 2 ))C^(a s (Q 2 )) (31) 

there is no place for the singlet-type contributions to the QCD generalization of the Crewther 
relation for the 3-point function of the axial-vector-vector non-singlet currents. The men- 
tioned in the item 4 cancellation of the light-by-light-type order a^-diagarms, which produce 
the colour structure (d abc ) 2 , is nothing more than the consequence of Eq.(31). The same 
equation is insuring this product from the contribution of the 3-loop triangle-type singlet 
diagrams, calculated in Ref. [47] and re-calculated with the same result in the works of Ref. [48] 
(for the details see Ref. [49]). 



n 



Three last "wonders" are leading to the striking indications on the factorization of the 
term (f3^(a s )/a s ) in the NNLO generalization of the Crewther relation of Eq.(28). It should 
be stressed, that we prefer to think about the theoretical origine of the factorization of this 
special term, but not simply of the 2-loop approximation of the QCD /3-function p( 2 \a s ), since 
in case of the cancellation of of a s in the denominator of the conformal symmetry breaking 
term {(5^ 2 \a s ) / a s ) by the extra power of a s in the square bracket of the r.h.s. of Eq.(28), we 
will get in it free Casimir operator CV, which is not multiplied by the coupling constant a s . 
In its turn, this will contradict the general property of the gauge structure of the concrete 
expressions for the perturbative contributions, which presumes the appearence of the Casimir 
operators only in the combinations, proportional to Cpa S) Cpa 2 , CfCacl 2 , CpTfNfa 2 ., etc. 

Another yet non-proved property of the QCD generalization of the Crewther relation is 
following from the "wonder" N6. Indeed, the absence of the NNLO contribution to A n _ c of 
the terms, proportional to fi^al is giving us the idea, that the explicite expression of the A n _ c - 
factor can not contain the contribution, proportional to ({3^ (a s ) / a s ) 2 . Though we cannot yet 
be sure, that the factor (f3(a s )/a s ) can be factored out of A n _ c beyond the NNLO, it seems 
to the authors of Ref. [5] most likely, that at any given order af {K > 1), one will encounter 
in A n _ c only the coefficients f3i\i < K, multiplied by linear combinations of colour factors. 

In fact we think, that the hypothesis A n _ c ~ ((3(a s )/a s ) merits close attention. In Ref. [6] 
the attempt to study the theoretical consequencies of this hypothesis in more detail was made. 
Using the assumption, that the factorized factor, which is appearing in the r.h.s. of Eq.(28) 
has the origin, similar to the one of the measure of the conformal invariance breaking within 
perturbation theory framework ((3(a s )/a s ) in the explicite expression for the anomaly of the 
energy momentum tensor it was proposed in Ref. [6] to rewrite Eqs.(16) for the QCD 
expressions of the tensor structures of the 3-point function of Eq.(2) as 

Uq\p 2 ) = (K(a s ) + ^v 2 (q 2 ,p 2 ,a s )y 2 ~ l ^( q 2 ,p 2 ) , (32) 

where v 2 and V3 are dimensionless functions, which can be constrained from the identity 

q 2 ^-M<l 2 ^ 2 ) = -P 2 4-2^ 2 ,P 2 )-U<l 2 ,P 2 ) • (33) 
aq z aq z 

This equation can be obtained from the derived in Ref. [23] Ward identity 

-6(<?V) =q%(q\p 2 )+p 2 Uq\p 2 ) (34) 

after noting that according to the Adler-Bardeen theorem the function £,i{q 2 ,p 2 ) is simply the 
unrenormalizable number. Keeping in mind Eqs.(32), which were written down after using the 
still non-proved assumptions that the factor (/3(a s )/a s ) is (a) indeed factorized in the QCD 
generalization of the Crewther relation and (b) it can really manifest itself in the perturbative 
expressions of £,2(q 2 ,p 2 ) and £,3(q 2 ,p 2 ) in the form, suggested in Eqs.(32) it is possible to show, 
that in QCD the Crewther product of Eq.(17) takes the following form ||: 

C B]P MQ 2 ))Cd S MQ 2 )) = 1 + —r{a s ) (35) 

where r(a s ) is a polynomial in powers of a s , which is not fixed in the approach of Ref. [6]. 
Therefore, the non-proved assumptions (a) and (b) are closely related. 
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Of course, a lot of work should be still done in order to find the really proved theoretical 
support in favour of these assumptions. Indeed, it is still necessary to understand on the di- 
agrammatic language the origine of the appearence of the term, proportional to {(3^ 2 \a a ) / a s ) 
in the discovered in Ref.[5] QCD generalization of the Crewther relation (see Eq.(28)). In 
view of the presented in Ref.[6] considerations we think, that the analysis of this problem 
should be started from the explicite calculations of the NLO perturbative QCD corrections 
to the invariant amplitude ^2(<? 2 ,P 2 ) of the 3-point function of Eq.(l). The second non-solved 
question is related to the necessity of the understanding, whether the factor (/3(a s )/a s ), pre- 
sumably related to the property of the conformal symmetry breaking in massless QCD, is 
indeed factorized in the generalized Crewther relation in all orders of perturbation theory. We 
believe, that the attraction of the formalizm of the non-conformal Ward identities, developed 
in Ref . [50] , can be rather useful for the study of this interesting problem. 

4 The generalzied Crewther relation and the "commen- 
surate scale relations" 

It should be stressed that the QCD generalization of the Crewther relation of Eq.(28) was 
discovered in Ref. [5] using the NNLO MS-scheme results for C% s (a s ) ]I§ [LJ and the similar 
ones for CBj P {a s ) BUI- However, it is known that on the contrary to the QED on-shell scheme, 
which is distinguished by the kinematics of the experimental measurements of the magnetic 
moments of electron and muon, the MS'-scheme is, regorously speaking, not physical one 
and is chosen as the reference scheme in the QCD phenomenology only by the convention 
between theoreticians and experimentalists. In view of this the number of approaches of 
dealing with the existing problem of fixing the scheme-dependence uncertainties in QCD was 



proposed. Among these methods are the principle of minimal sensitivety [^T|], the effective 
charges approach ]52"| , [5B| (which is known to be a posteriori equivalent to the scheme-invariant 
perturbation theory, developed in Refs. [54,55]) and the BLM approach [|56|j , generalized to the 
NNLO level in Ref. [57] (see also Ref. [58]). 

Following the work of Ref. [21] we will show that it is possible to apply the methods of 
Refs. [51-56] in order to rewrite the generalized Crewther relation of Eq.(28) in the different 
form and to derive the variants of the obtained in Ref. [59] so-called "commensurate scale 
relations". Let us first define the effective charges of the coefficient functions C B s (a s ) and 
C B j P (a s ) as 

C% s (a s ) = l + MQVAf,) , 

where an = ^^af and ctsjp = ^^af™ can be related to the MS'-scheme results using the 
following equation 

af^ D \Q 2 /A 2 BMD) ) = a s (Q 2 /Ak) + ( A m + B m (3»)a]{Q 2 / A^) 

+ {C m + D m (3 + E m f3l + Br^aliQ 2 /^) . (37) 

Here A 2 B - p ^ = Aj^exp ( Al(2) + p Q 1{2) ^° ^ are the effective scales of CBj P (a s ) and C^ s (a s )-coeflicient 

functions, (3$ and (3\ are the scheme-invariant coefficients of the QCD /3-function (see Eq.(20)) 
and the exact expressions of the terms Ai, Bi,C{, Di, Ei (i = 1,2) can be extracted from the 
NNLO results of Eqs.(25),(24). Two effective charges af^ p and af can be related as 
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af JP (Q 2 /A| JP ) = of (QV^n) + (An + B 12 [3 )(a?(Q 2 /A 2 D ) 2 

+ (C12 + D l2 (3 + E 12 [3 2 + B 1 M{a°(Q 2 /A 2 D )) 3 (38) 

where A 12 = A 1 -A 2 , B 12 = B 1 -B 2 , C l2 = C l -C 2 -2(A l -A 2 )A 2 , D 12 = D 1 -D 2 -2{A 1 B 2 + 
A 2 Bi) + AA 2 Bx, Ei 2 = Ei — E 2 — 2(B\ — B 2 )B 2 are the scheme-independent coefficients. 

The basic aim of the applications of the BLM criterion ]56]| to the series of Eq.(38) is 
the construction of the iV/-independent perturbative scheme-independent expansion, which 
can play the role, similar to the one of the Baker- Johnson scheme- independent Fi-funcion in 
QED. At the NNLO level this aim can be achieved with the help of the developed in Ref.[57] 
single-scale generalization of the BLM method, which leads to the following redifmition of the 
scale of the af -effective coupling constant from Q 2 to Q* ^jj : 

ln(Q* 2 /Q 2 ^J = -B 12 + [(3 (B 2 2 - E l2 ) + 2A 12 B 12 - D 12 ]af (Q* 2 / k 2 D ) . 

It should be stressed, that this procedure is effectively eliminating the dependence of the r.h.s. 
of Eq.(38) from the /^-coefficient. Moreover, since the coefficient before the /?i-term in the 
NNLO contribution to Eq.(38) can be chosen to be equal to the coefficient B\ 2 of the NLO 
correction, the absorption of the proportional to Nf (and thus to (3q) NLO term into the scale 
Q automatically leads to the nullification of the proportional to f3\ NNLO contribution in 
this relations bewteen effective charges. 

Taking now into account the concrete expressions for the effective charges af and af^ p (see 

Eqs.(24),(25)) one can find [|TJ ln(Q* 2 /Q 2 ) = |-4C(3) + af (Q*) + f C(3) - 16C 2 (3))/? + 



§C A - §CUC(3) - 1 -§C F - f C(3)C F + 20((5)C F 



and 



a Bjp {Q) = a D {Q)-aj 3 (Q)+a 6 D {Q) + - • (39) 

Equation (39) is representing the example of the single-scale variant of the "commensurate 
scale relation" of Ref. [59]. It is the consequance of the following variant of the generalized 
Crewther relation of Eq.(28): 

C Bjp (a^(Q 2 ))C% s (a?(Q* 2 ) = 1 (40) 

where the non-conformal term A n _ c in Eq.(28) is absorbed into the scale of the coefficient func- 
tion Cp S . It should be stressed, that the property of the factorization of the term {(3^ 2 \a s ) / a s ) 
in the expression of Eq.(28) for A n _ c , discovered in Ref. [5], turns out to be the necessary 
and the sufficient condition, which allowed the authors of Ref. [21] to rewrite the generalized 
Crewther relation of Eq.(28) in the form of the geometric progression of Eq.(39). Note also, 
that the less convenient for practical applications multi-scale variant of the "commensurate 
scale relation" of Eq.(39), which was previously derived in Ref. [59], is also the consequence of 
the definite variant of the generalized Crewther relation of Eq.(28). 



5 The generalized Crewther relation and the experi- 
ment 

After reading the previous Sections one can be interested in getting the understanding whether 
it is possible to build the bridge between the pure theoretical studies presented above and the 
existing experimental data for the characteristics of the e + e~-annihilation and deep-inelastic 
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processes, which enter in the Crewther relation and its QCD generalizations. The first at- 
tempts to analyze this problem was made in Ref. [21] and in the closely related work of Ref. [58] . 
In this Section we will follow the original considerations of Ref. [21], leaving the comments to 
the outcomes of Ref. [58] for the possible future presentation. 

Let us first mention, that while the deep-inelastic scattering sum rules are measured in the 
Eucledian region, the real experimental information for the basic chracteristic of the e + e~- 
annihilation channel R e + e -(s) = cr to4 (e + e~ — > hadrons)/a(e + e~ — > is coming from 

the measurements in the Minkowskian region. If one is sufficiently far from the resonance 
production thresholds, it is possible to relate the perturbative expression for R e + e ~(s) with 
the coefficient function C^ M (Q 2 ) by the following relation 



R e +e~ 



3E/Q/ f' s+i 'dr 



2ni 



-C 



EM i 

d \ 



(r)) 



(41) 



In general this procedure results in the appearence of the 7r 2 -like terms in the coefficient 
function of the i?-ratio starting from the NNLO-level (see e.g. Ref[60]). In fact they can be 
taken into account in the formulae of Eqs.(37),(38) using the following shifts E 2 — > E 2 — vr 2 /3, 
E 12 — > Ei 2 + 7r 2 /3. The higher-order ^-contributions to R e + e -(s) were calculated explicitely 
in Refs. [61,62]. However, since at the present stage we are interested in the consequencies of 
the NNLO generalizations of the Crewther relation, we will not sum-up these higher-order 
^-contributions (like it was proposed to do e.g. in Ref. [63]), but trancate the corresponding 
perturbative series in the time-like region at the NNLO level. In order to obtain the relation 
between the time-like scale \/s* of the _R e + e --ratio and the space-like scales of the deep-inelastic 
scattering sum rules it is not only necessary to take into account the NNLO ^-contributions 
to the i?-ratio, but to perform the replacements of the effective charges af (Q*) — > (\/s*) as 
well. Than the resulting relation reads [[2"T|]: 



ln{^) =-I+4C(3)-af(0 



I + f CO) - 16C 2 (3) - f>o 



^C A - ^C A C(3) - ^C F - yC(3)CV + 20C(5)CV 



(42) 



The variants of the generalized Crewther relation can be then written down in the following 



form 5T 



1 



C 



Bjp{ 



(Q 2 ))R e 



SI(„R, 



O) 



(43) 



1 



C G Ls(af LS m)R e+e -(a? (s*)) = 1 + C? (af {s*)) + C^ LS (a^(Q 2 )) 



SIf„Ri 



iSI 



MLS/ 



3E f Qf 



where Cf 1 abd Cq LS are the singlet-type contributions to C 



EM 1 
D \ 



and CcLsifls) which are 



defined in Eq.(23) and Eq.(27) respectively. In fact, since the numerical values of these con- 
tributions are very small, it is reasonable to neglect them in the phenomenologically oriented 
discussions and thus assume, that C^ M (a s ) ~ C^ s (a s ) and CBj P {a s ) ~ Cgls{&s)- 

It should be stressed, that the theoretical expressions of Eq.(43) are derived in the frame- 
work of the perturbation theory and do not involve the non-perturbative contributions to 



C 



EM 1 
D \ 



(and thus i?-ratio), and to the deep- inelastic scattering sum rules, theoreti- 



cally calculated in Ref. [65] and numerically estimated using the QCD sum-rules formalizm 
in Ref. [66]. In fact it is known, that these contributions are very important for the analy- 
sis of the low-energy experimental data for the i?-ratio p7| and the Gross-Llewellyn Smith 
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and the polarized Bjorken sum rules (see Refs. [68,69] and Refs. [70,71] respectively). There- 
fore, in the estimates, which have the aim to clarify the experimental status of the per- 
turbative QCD generalizations of the Crewther relation, it is necessary to chose the scales 
in the regions of energies, where the non-perturbative effects can be safely neglected. For 
the e + e~-annihilation i?-ratio the lowest energy region, which is satisfying this criterion, is 
4 GeV < y/s < 8 GeV. It is typical to the case of taking Nf = 4 numbers af active flavours 
and is lying above the thresholds of the production of the cc-bound states y/s ~ 3 GeV and 
beyond the thresholds of the production of 66-bound states y/s ~ 10 GeV. The detailed QCD 
analysis of the e + e~-data for the i?-ratio in this region was made previously in the number 
of the works on the subject (see e.g. Refs. [72,73]). In Ref.[21] the result of the most recent 
similar fit of Ref. [74] was used, which provide the authors of Ref.[21] with the constarint 
(1/3 £/ Q 2 f )R e+e - (y/s* = 5.0 GeV) « 1.08 ± 0.03 and thus af (V^ = 5.0 GeV) « 0.08 ± 0.03. 
Equations (42), (43) then imply that the corresponding estimates for the coefficient functions 



of the sum rules lie on the segment line which connects the following three values [21 



C BjP (Q = 11.13 GeV) « C GLS {Q = 11.13 GeV) « 0.952 , 

CbM = 12.33 GeV) « C GLS {Q = 12.33 GeV) « 0.926 , (44) 

CbjpIq = 13.53 GeV) « C GLS {Q = 13.53 GeV) « 0.900 . 

The corresponding values of the effective coupling constants lie in the line segment, which 
connects the three related points 

af jp (Q = 11.13 GeV) « af LS (Q = 11.13 GeV) » 0.048 , 

af>(Q = 12.33 GeV) « af LS (Q = 12.33 GeU) w 0.074 , (45) 

af ip (Q = 13.53 GeV) « af L5 (Q = 13.53 GeV) « 0.1 . 

The appearence of three corresponding values of Eqs.(44) and Eqs.(45) are related to the 
uncertainties in the definition of af(y/s = 5 GeV), which are translated by Eq.(42) into three 
related numbers. The predictions of Eq.(45) could be tested experimentally. 

It should be stressed, that at present the measurements of the polarized Bjorken sum 
rule are allowing to obtain its precise experimental values in the regions of energies of over 
3 GeV 2 (see Ref. [75]) and 10 GeV 2 (see Ref. [76]). The recent measurements of the Gross- 
Llewellyn Smith sum rule are also preformed at relatively small values of Q 2 (see Ref. [69]). 
However, for the estimates, aimed to the comparison with the results of Eqs.(45), it was 
proposed in Ref. [21] to use the results of the theoretical extrapolation [|77j of the available 
experimental data of the CCFR collaboration |78| to the wide region of energies. The indirect 
determination of the values of the Gross-Llewellyn Smith sum rule for 3 GeV 2 < Q 2 < 
500 GeV 2 of Ref. [77] gives the value 

a ? LS (Q — 12.25 GeV) w 0.093 ± 0.042 . (46) 

This interval crosses the line of the estimates of Eq.(25). To the point of view of the authors of 
Ref. [21] this fact gives empirical support for the generalized Crewther relation, written down 
in the form of Eq.(40). 

In order to clarify the meaning of these numbers it seems to us rather instructive to return 
to the commonly used language of the MS'-scheme. Let us first analyse the question of the 
extraction of the value of Aj^g from the fits of the e + e~-data of Ref. [74]. The results of the 
NNLO calculations of Refs. [10,11] relate the effective charge of the e + e~-annihilation to the 
coupling constant a s in the MS'-scheme [1 + 1.524a s - 11.52e#. Here N f =A is taken 

and the small contribution of Gf 7 (a s ) is neglected. Using this relation, we find that the used 



by us experimentally motivated number for af(y^s = 5 GeV) corresponds to the following 
value of the coupling constant a s in the MS'-scheme: a s {yfs = 5 GeV) = 0.238+o;og 2 . It is 
known that the coupling constant a s can be expressed through the QCD scale parameter 
as 

1 AlnlntQVAf^) 



4vr (3 \n(QyAy) /3 2 ln 2 (QVA|^) 
g In 2 ln(g 2 /A|^) - (3 2 In hi(Q 2 /A^) + frfa - (3 2 

Combining Eq.(47) with the obtained above value of a s {\fs) (normalized to Nf = 4 num- 
bers of active flavours) we find the following interval of the values of the parameter A|^ = 
410i33Q MeV. Evolving this result through the threshold of the production of the 66-bound 
M = 2m b ps 9 GeV using the approximate formula of Ref. [79] we get A^ = 275^0 MeV 
and 

a s {M z ) = 0.123l°;°27 • (48) 
This value is compatible with the result a s (Mz) = 0.124 ±0.021, which comes from the recent 



detailed fit of the available data in e + e~-annihilation from yfs = 20 GeV tp 65 GeV |80 
(for the detailed discussions of this result see the recent review of Ref. [81]. Thus the current 
measurements of R e + e - suffer from sizable experimental uncertainties. 

The similar situation also holds in the case of the existing experimentally motivated data 
for the Gross-Llewellyn Smith sum rule at the energies, higher than Q 2 > 5 GeV 2 (see, for 
example, the results of the recent analysis of Ref. [69]). In order to demonstrate explicitely 
that the used outcomes of the extrapolation-extraction of the Gross-Llewellyn Smith sum rule 
value at Q 2 = 150 GeV 2 (see Ref. [77]) has definite theoretical and experimental uncertainties, 
we transform the results of Eq.(46) into the MS'-scheme and obtain a s (Q = 12.25 GeV) ~ 
0.220lo;ogg. This estimate might be weekly senistive to the change of Nf from Nf = 4 (which 
is typical to the fits of the deep-inelastic data) to Nf = 5, which is more appropriate to the 
scale Q 2 = 150 GeV 2 . Therefore, we will extract from the high-energy results of Ref. [77] the 
value of Ai±. and will get 

a s (M z ) = 0.119™ ■ (49) 

One can see that this estimate has larger uncertainties, than the value a s (Mz) = 0.109 ± 
0.003(stat) ± 0.005(syst) ± 0.003(thor), recently extracted in Ref. [82] from the analysis of 
the CCFR data for xF 3 structure function using the information about the NNLO correc- 



tions to the coefficient functions [Q and the results of the recent calculations of the NNLO 
corrections to the anomalous dimensions of the non-singlet operators |j84| . Therefore, the cur- 
rent experimentally-motivated estimate of the value of the Gross-Llewellyn Smith sum rule 
at Q 2 = 150 GeV 2 (see Ref. [77]) also suffers from the definite theoretical and experimental 
uncertainties. 

These discussions are giving the additional arguments in favour of the physical conclusions 
of Ref. [21]: in order to check the consequencies of the perturbative generalizations of the 
Crewther relation at more high confidence level it can be rather helpul, first, to reduce the 
experimental error of the measurements of R e + e - at ^fs ps 5 GeV and, secondly, to have more 
precise information on the value of the Gross-Llewellyn Smith sum rule or of the Bjorken 
polarized sum rule at Q 2 ~ 150 GeV 2 . The first problem can be attacked after starting the 
operation of the c — r-factory, while the possible future study of the deep-inelastic scattering 
with both polarized electron and proton beams can open the window for direct measurements 



of the polarized Bjorken sum rule at high momentum transfer [85 
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Another interesting, to our point of view, proposal is to try to measure the value of the po- 
larized Bjorken sum rule (or the Gross-Llewellyn Smith sum rules) at the scale Q 2 = 25 GeV 2 . 
It can be useful for the study of the status of the perturbative generalization of the Crewther 
relation, written down in the form of Eq.(28), derived in Ref. [5]. We hope to return to the 
more detailed analysis of this problem in our possible future work. 

6 Is there any generalzation of the Crewther relation 
for the 3-point function of singlet axial- non- singlet 
vector-vector currents? 

As was already explained in Sec. 2, the classical Crewther relation was originally derived in 
Ref. [3] for the 3-point function of the axial- vector- vector non-singlet currents. For the sake 
of completeness let us now discuss the case of the analogouse 3-point function with the axial 
singlet current : 

T;Up> ?)=»*/< 0\TAMV2{x)V*{0)\0 > e***«dxdy (50) 

where = i^j^jsip. In this case it is also possible to try to repeat the considerations of Sec. 2 
and to think about the possibility of the derivation of the Crewther-type relations. Indeed, 
keeping the singlet structure in the operator-product expansion of the two non-singlet vector 
currents, one can get 

% J TV^e ipx dx\ p 2^ « C^' p ab A p (o) + other structures . (51) 

In analogy with Eq.(10) one finds that 

~ i«F%^Cg(a.) . (52) 

The coefficient function Cjjj(a s ) is entering in the definition of the so-called Ellis- Jaffe sum 
rule as |56 

EJ(Q 2 ) = f 1 g p 1 (n \x,Q 2 )=(±±\g A \ + ±a 8 )c% J s (a s ) 



(53) 

+ |AS(/i 2 ) x AD x C|j(a s ) 

where Cff (a s ) = C Bjp (a s ), AD = expff^^ 'Zj^dx) and \g A \ = An- Ad, a 8 = Au + Ad- 

2As, AS = An + Ad + As, Am, Ad, As can be interpreted as the measure of the polarization 
of quarks in a nucleon, C^j(a s ) = C^ p (a s ) and 7 S7 (a s ) is the anomalous dimension of the 
singlet axial current, which can be calculated within the MS'-scheme provided the additional 
finite renormalization of the singlet axial current is made: 

An -> Zl l Ap = A„ (54) 

It was shown in Ref. [38] that the finite renormalization constant Z^ 1 is different from the 
defined in Eq.(21) "non-singlet axial charge" due to the additional contribution to Z$ of 
the 3-loop light-by-light-type triangle diagrams. At the 2-loop order this difference is Z^ 1 = 
Z^ s + -^CpNfa 2 [^1 - The similar difference is also appearing in the coefficient function 



CfJ(a s ) for the Adler .D-function of two singlet currents, which can be constructed from the 
following correlation function 



i 



J <0\TA^x)A v (0)\0>e i * x dx = IL s J,tf) . (55) 



Staring from the 3- loop order the coefficient function C^/(a s ) differs from the coefficient func- 
tion Cp (a s ) by the light-by- light-type diagrams, first calculated in Ref. [47] in the limit of the 
heavy top-quark mass m t : 



C b D 1 (a s )=C% b (a s )+AC b D I (a s ) . (56) 

In the case of the lighter quarks m q « m t the result for the 3-loop light-by-light-type con- 
tribution can be extracted from the first paper of Ref. [49]. 

Let us now return to the discussion of the singlet contribution to the Ellis- Jaffe sum rule. 
The NLO corrections to the coefficient function C^j(a s ) was calculated in Ref. [86]. It is 
different from the result of the calculation of Ref. [12] due to the appearence of the light-by- 
light-type singlet diagrams 

CfjK) = C™(a B ) + AC|S(a s ) = C Bjp {a s ) + ACf5(a.) (57) 



where at the NLO level AC^j(a s ) = (C(3) + ±C F T f N f )a 2 s This result is in agreement 
with the calculation of Ref. [83]. The anomalous dimension 7 S7 (a s ) was calculated in the 
MS-scheme at the 3-loop order in Ref. [38]. In Ref. [27] the order a^-term was related to 
the calculated in Ref. [26] light-by- light-type diagrams, contributing to the 3-point function of 
Eq.(50). 

Let us now consider the conformal-invariant limit of all discussed in this Section results. 
In this case we have the Crewther-type relation 

C^MQ 2 ))^ 1 (a s (Q 2 ))| c „,,= 1 . (58) 

We do not still know, whether there exist its QCD generalization, analogous to the one of 
Eq.(28), discovered in Ref. [5]. From the general grounds we expect, that both sides of Eq.(58) 
should be multiplied by the anomalous dimension term AD, defined in Eq.(53). It is interesting 
to study the problem of the possibility of the existence of any relations between the extra 
light-by-light-type terms ACfj 7 (a s ), ACjj(a s ) and the light-by-light-type QCD contributions 
to the r.h.s. of Eq.(58). The most surprizing fact is that the explicite analytical expression for 
ACj j(a s ) contains the £(3)-term P6[ , while the 3-loop correction ACf/(a s ) does not contain 
it (see the first work of Ref. [49]). We think that the calculations of the NNLO corrections 
to the Ellis- Jaffe sum rule [|8^] might be useful not only for the comparison of their results 
with the scheme- invariant estimates of Ref. [88] , but for the analysis of the possibility of the 
existence of new "wonders" in the non-conformal deviations from the Crewther-type relation 
of Eq.(58). 
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